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We study the spontaneous emission, the absorption and dispersion properties of a A-type atom 
where one transition interacts near resonantly with a double-band photonic crystal. Assuming an 
isotropic dispersion relation near the band edges, we show that two distinct coherent phenomena 
can occur. First, the spontaneous emission spectrum of the adjacent free space transition obtains 
'dark lines' (zeroes in the spectrum). Second, the atom can become transparent to a probe laser 
field coupling to the adjacent free space transition. 

(3 , I. INTRODUCTION 

O; 

^^ , It has been well established over the years that spontaneous emission depends not only on the properties of the 

CIh" excited atomic system but also on nature of the surrounding environment and more specifically on the density of 

1^ , electromagnetic vacuum modes. Purcell |ll was the first to predict that the rate of spontaneous emission for atomic 

transitions resonant with cavity frequencies could be enhanced due to the strong frequency dependence of the density of 

\^ _ modes. Employing a similar idea, Kleppner B] showed that spontaneous emission in a waveguide may be suppressed 

CN ■ below the free space level due to the singular behaviour of the density of modes near the fundamental threshold 

frequency of the waveguide. 

The above studies attracted recent attention when it was realized that periodic dielectric structures could be 
engineered in such way that gaps in the allowed photon frequencies may appear p]. Since then, photonic band gap 
f^ , (PBG) materials exhibiting photon localization have been fabricated initially at microwave frequencies and more 
^— I ■ recently at near infrared and optical frequencies. The study of quantum and nonlinear optical phenomena in atoms 
^^ ' (impurities) embedded in such structures, lead to the prediction of many interesting effects 0]. As examples we 
^D mention the localization of light and the formation of single-photon ^-0] and two-photon |^ 'photon-atom bound 
>^ ■ states', suppression and even complete cancellation of spontaneous emission &-n2], population trapping in two-atom 
^^ : systems p2(| , phase dependent behaviour of the population dynamics ||l^,|lj| , enhancement of spontaneous emission 
f| ' interference |15|"p7|, modified reservoir induced transparency |18|-p0| and transient lasing without inversion [ pi[ , 
Qh. occurrence of dark lines in spontaneous emission ||2^,^ and other phenomena [p4|-p6|. In addition there is also 
.«_> ■ current interest with regard to the feasibility of the observation of either the quantum Zeno effect p^ or the quantum 
H Anti-Zeno effect |g^ in modified reservoirs, such as the PBG. 

In this paper we study the spontaneous emission and the probe absorption/dispersion spectrum in a A-type sys- 
tem, similar to the one used in previous studies ||l^;|l^,|l^,^ 22 2jj22l , with one of the atomic transitions decaying 
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spontaneously in a double-band [g/7| PBG reservoir. We show that the spontaneous emission spectrum can exhibit 

.'fl ' 'dark lines' (zeroes in the spectrum at certain values of the emitted photon frequency) in the spontaneous emission 

K^ , spectrum in the free space transition. In addition, the atom becomes transparent to a probe laser field coupling to 

^ ' the free space transition, which exhibits ultra-slow group velocities near the transparency window. The dependence 

^ ' of these phenomena on the width of the gap is examined and a comparison with the single-band case results is made 

This article is organized as follows. In section II we apply the time-dependent Schrodinger equation to describe 
the interaction of our system with the modified vacuum, calculate the spontaneous emission spectrum in the free 
space reservoir and discuss its properties along with the differences from the single-band case |^,0. In section III 
using the time-dependent Schrodinger equation we calculate analytically the steady state linear susceptibility of the 
system. Results of the absorption and dispersion of a probe laser field in this system and their comparison with the 
single-band case are also presented in section III. Finally, we summarize our findings in section IV. 

II. EQUATIONS AND RESULTS FOR THE SPONTANEOUS EMISSION SPECTRUM 

We begin with the study of the A-type scheme, shown in Fig. 1(a). This system is similar to that used in previous 
studies |10|jl^Jlq , pi| . p2| , p4] . ^ . The atom is assumed to be initially in state |2). The transition |2) ^^ |1) is taken to 
be near resonant with a modified reservoir (this will be later referred to as the non-Markovian reservoir), while the 
transition |2) <-^ |0) is assumed to be occurring in free space (this will be later referred to as the Markovian reservoir). 
The spectrum of this latter transition is of central interest in this section. The Hamiltonian which describes the 



dynamics of this system, in the interaction picture and the rotating wave approximation (RWA) , is given by (we use 
units such that ?i = 1), 

H^Y. 5Ae"*('^^-"-)*|2) {0\ax + J2 .9.e-^("''-'^-)*|2) (l|a, + H.c. . (1) 

A K 

Here, g^ denotes the couphng of the atom with the modified vacuum modes (k) and g\ denotes the couphng of the 
atom with the free space vacuum modes (A). Both couphng strengths are taken to be real. The energy separations of 
the states are denoted by Wy — uji — loj and uj,^ {^x) is the energy of the k (A)-th reservoir mode. 

The description of the system is done via a probabihty amphtude approach. We proceed by expanding the wave 
function of the system, at a specific time t, in terms of the 'bare' state vectors such that 

|^(i)) = 62 W|2, {0}) + Y. 6.(i)|0, {A}) + Y ^'^Wll' {4> • (2) 

A K. 

Substituting Eqs. (|l|) and ^ into the time-dependent Schrodinger equation we obtain 

Mt) - Y. 9xbx{t)e-'^^^-^'°^' + Y 9.b.{t)e~'^^'^-^''^' , (3) 

A K 

ibxit) = gxb2it)e'^^'-'^'°'>\ (4) 

i&.(i)=ffK&2(i)e*^"'= -"-'*. (5) 

We proceed by performing a formal time integration of Eqs. (Q) and (g) and substitute the result into Eq. (Q) to 
obtain the integro-diffcrcntial equation 



hit) = -fdt'hit')Y9le-'^^'-^-'''^^'~''^-fdt%it')Y9y 
Ja , JO 



i(w„-cJ2l)(t-t') 



(6) 



Because the reservoir with modes A is assumed to be Markovian, we can apply the usual Weisskopf-Wigner result 0] 
and obtain 

Y (7^e-''('^^-'^™)(*-*') = ^8{t - t') . (7) 

A 

Note that the principal value term associated with the Lamb shift which should accompany the decay rate term has 
been omitted in Eq. ffl). This does not affect our results, as we can assume that the Lamb shift is incorporated 
into the definition of our state energies. For the second summation in Eq. (^) , the one associated with the modified 
reservoir modes, the above result is not applicable as the density of modes of this reservoir is assumed to vary much 
quicker than that of free space. To tackle this problem, we define the following kernel 

with /3 being the atom- modified reservoir resonant coupling constant. The above kernel is calculated using the 
appropriate density of modes p{lo) of the modified reservoir. Using Eqs. (0) and (g) into Eq. (||) we obtain 

nt 

Ht) = -lb2{t) - J dt'b2{t')K{t^t'). (9) 

The long time spontaneous emission spectrum in the Markovian reservoir is given by S{5x) oc \b\{t — * 00) p, with 
5\ = uj\ — W20- We calculate b\(t ^ 00) with the use of the Laplace transform of the equations of motion. Using Eq. 
(0) and the final value theorem we obtain the spontaneous emission spectrum as S{5\) ex 7I Mvas^-iS), -B2(s)P, where 
-82(5) — j^ dte^^*b2{t) is the Laplace transform of the atomic amplitudes 62(0 ^^'^ ^ is the Laplace variable. This 
in turn, with the help of Eq. (S), reduces to 

SiSx) cc ^^ 2 . (10) 

-iSx + 7/2 + K{s -^ -iSx) 

Here, K{s) — /„°° dte'^'^^Kit) is the Laplace transform of K{t), which yields from Eq. (g), 



K{s) = /33/2 f diU 



S + i{LJ ~ UJ21) 



(11) 



Therefore, in order to calculate the spontaneous emission spectrum in the Markovian reservoir we need to calculate 

k{s). 

We consider the case of a double-band isotropic model of the photonic crystal, which has an upper band, a lower 
band and a forbidden gap, the near the edges dispersion relation p7| is approximated by 



with An w ujgn/KQ, (ji = 1, 2). Then the density of modes reads, 



(12) 
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y/LU - LUg2 



Q{UJ - UJg2) 



(13) 



where LOgn is the gap frequency and O is the Heaviside step function (as shown in Fig. 1(c)). Then, from Eq. ( [ll[ ) we 
obtain 



k{s) 



I3^'\ 



yjs + i6. 



ffl 



^3/2 



i5, 



s2. 



(14) 



with 5gn = LOgn ~ ^21. 

We use the formulae obtained above and calculate the spontaneous emission for several parameters of the system. 
From Eqs. (O) and ([l4h is obvious that the spectrum exhibits two zeros (i.e. predicts the existence of two dark 
lines), at 5\ — 6gi and Oa = <^g2- This is purely an effect of the density of modes of Eq. (p^), and the non-Markovian 
character of the reservoir. In the case of a Markovian reservoir the spontaneous emission spectrum would obtain the 
well-known Lorentzian profile and no dark line would appear in the spectrum. The behaviour of the spectrum is 
shown in Fig. 2. In Fig. 2(a) we show an asymmetric case where the right hand side (rhs) threshold detuning {Sg2) is 
kept constant while the left hand side (Ihs) threshold detuning (Sgi) is changed. There are three peaks and two zeroes 
in the spectrum, one pronounced peak in the center and two lower peaks in the sides. A s Sqi g oes further to the sides 
the Ihs peak is suppressed and the spectrum approximates that of the single band case [ |10[|ll| ] which is shown in Fig. 
3. That is because the Ihs singularity in the density of the reservoir modes, Fig 1. is far detuned from the atomic 
transition. 

In Fig. 2(b) we show a symmetric spectrum, which is the case when the atom is placed in the middle of the gap. 
This spectrum also has two zeroes and three peaks. In this case when the threshold detunings increase (i.e, the width 
of the gap increases) the spectrum approximates the Lorentzian shape. This is a result of the complete cancelation 
of any emission in the modified reservoir as there are no modes available for a large range of frequencies around the 
relevant atomic transition. 



III. EQUATIONS AND RESULTS FOR THE PROBE ABSORPTION/DISPERSION SPECTRUM 

The aim in this section is to investigate the absorption and dispersion properties of our system for a weak probe 
laser field. As the probe laser field is assumed to be weak the probability amplitude approach can be employed 
for the description of the system. The Hamiltonian of the system, in the interaction picture and the rotating wave 
approximation, is given by. 



i7 = fie""|0)(2| 






gxe 



-i{uJx~^2o)t 



|2)(0|aA 



,.9^6 



— i{u)K—u)2l)t 



|2)(1| 



H.c. 



(15) 



Here, Q, is the Rabi frequency (assumed real for simplicity) and 5 = uj — u>2o is the laser detuning from resonance with 
the |0) <-^- 1 2) transition with cu being the angular frequency of the probe laser field. We note that we are interested 
in the perturbative behaviour of the system to the probe laser pulse, therefore we can eliminate the Markovian 
decay modes using the method presented in the previous section and study the system using the following effective 
Hamiltonian 



H 



Qe'^'\0){2\ + J2 



5k e 



— i(uj^—ui2i)t 



|2)(l|a« + H.c. 



^lmm 



(16) 



The wavefunction of the system, at a specific time t, can be expanded in terms of the 'bare' eigenvectors such that 
m)) = co(t)|0, {0}) + C2(i)e-*'*|2, {0}) + J2 c« W|l' M) ' (17) 



and co(i = 0) = 1, C2(i = 0) = 0, CK(i = 0) = 0. Substituting Eqs. ( |16D and ( |17D into the time-dependent Schrodinger 
equation and ehminating the vacuum ampHtude c^it), we obtain the time evolution of the probabihty amphtudes as 



icoit) ^ nc2{t) , (18) 

iC2it) = Qcoit) - (6 + i^) C2(t) - I I dt'K'it - t')c2{t') , (19) 

with the kernel 

K\t -t')=Y^ ^2g-i(^„-a;2i-A-)(t-t') _ K{t)e'^' . (20) 

All the coupling constants {g^, (3, f2) are assumed to be real, for simplicity. 
The equation of motion for the electric field amplitude E{z, t) is given by [?], 

| + ^|)i^(.,t) = -.|-x(W^,^), (21) 

where x((5) is the linear susceptibility of the medium and Vg — c/[l + {Lo/2){dY{.e{x)ldLo)] is the group velocity of the 
laser pulse with the derivative of the real part of the susceptibility being evaluated at the carrier frequency. 

Since the transition |0) <-> |2) is treated as occurring in free space, the steady state linear susceptibility is given by 

X{5) = -^^^^co(i - oo)c^(t ^ oo) , (22) 

with A/" being the atomic density. 

Therefore, in order to determine the steady state absorption properties we have to solve Eqs. (Eq), (M) for long 
times. We assume that the laser-atom interaction is very weak {Q. <C /3,7) so that Co(i) ~ 1 for all times, and by 
using perturbation theory, Eqs. (|l8|) and ( pj| ) reduce to 

ic2{t) KVL-U + i^ C2{t) - i I dt'K'it - t')c2{t') . (23) 

We further assume that Q{z,t) is approximately constant in the medium and with the use of the Laplace transform 
we obtain from Eq. (|2^ ) 

C2{s) = -. — J . (24) 

s 6 + i-f/2 + iK{s) + is 



If 7 ^ then the terms inside the brackets of Eq. (24) have only complex, not purely imaginary, roots. Therefore we 



can easily obtain, using the final value theorem, the long time behaviour of the probability amplitude. 



C2{t ^ Qo) = lim[C2(s)] 



n 



5 + ?:7/2 + iK'{s -^ 0) 

= ^^ . (25) 

S + ij/2 + iK{-iS) 

Using Eqs. ( p^ ) and (|2^) we find that the steady state linear susceptibility of our system is given by 

X(^) = '"-^l^"!'' . (26) 



Therefore, the absorption and dispersion of the probe laser field are determined by the density of modes of the non- 
Markovian reservoir. In our case the linear susceptibility becomes zero at S = Sgi and S = Sg2, therefore the medium 
becomes transparent to the laser field. This result is in contrast with the case in which the transition |2) ^^ |1) occurs 
in free space, where the well-known Lorentzian absorption profile is obtained. Typical spectra are shown in Fig. 4, 
where only the Ihs threshold detuning is changed. Both the absorption and dispersion spectra are asymmetric and 
their shape depends critically on the dctunings. We note that Fig. 4(c) resembles as expected, the corresponding 
spectra in the case of a single band reservoir |Tq] which is shown in Fig. 5. 

We also present the case of symmetric spectra in Fig. 6, where the detunings are symmetrically changed i.e., with 
the transition frequency of the atom in the middle of the gap. Here the interesting effect is best shown in Figs. 6(a) 
and 6(b) where the transparency effect in the sides of the spectra is combined with the usual absorption/dispersion 
profile near the center of the spectrum. 

We note that besides the transparency effect which occurs for the predicted values of frequency, we note that the 
group velocity of the probe laser pulse is also effectively reduced near the transparency window due to the steepness 
of the dispersion curve. These results are related to that of electromagnetically induced transparency (FIT) in a 
four-level laser-driven tripod scheme [Q. However, KIT occurs through the application of external laser fields. Here, 
transparency is intrinsic to the system as it occurs due to the presence of a square root singularity of the density of 
modes at the band edges. 

IV. SUMMARY 

In this article we studied the spontaneous emission, absorption and dispersion spectra of a three-level atom embed- 
ded in a two-band isotropic photonic crystal. For the spontaneous emission we have shown that the spectrum exhibits 
two dark lines and can be significantly modified via the system parameters. The explicit dependance of the spectrum 
on the width of the gap and on various values of the atomic parameters was also analyzed. For the absorption and 
dispersion spectra, two transparency windows exist. Near the transparency windows a very slow group velocity for 
the laser pulse is obtained. In this case too, the spectra depend critically on the atomic parameters and the width of 
the gap. 
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FIG. 1. Figure (a) displays a three level, A- type atomic 
system. The thick dashed line denotes the coupling to the 
modified reservoir (PBG) and the thin dashed line denotes 
the background decay. Figure (b) shows the density of modes 
for the case of the single-band isotropic PBG model. Figure 
(c) shows the density of modes for the case of the double-band 
isotropic PBG model. 
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FIG. 2. The spontaneous emission spectrum S{S\) (in ar- 
bitrary units) for the double-band case and parameters 7 = 1 
and (a) (Sgi = —1, 5g2 = 0) (full curve), {Sgi = —2, 5g2 = 0) 
(dashed curve), [Sgi = —3, (5g2 = 0) (dot-dashed curve), (b) 
{Sgi = —1, Sg2 = 1) (fuU curve), {Sgi = —2, 552 = 2) (dashed 
curve), {Sgi — —3, (592 ~ 3) (dot-dashed curve). 
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FIG. 3. The spontaneous emission spectrum S{5\) (in ar- 
bitrary units) for the single-band case and parameters 7 = 1 
and Sg — (full curve); 5g = 1 (dashed curve); Sg — —1 
(dot-dashed curve). All parameters are in units of l3. 
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FIG. 4. The absorption and dispersion spectra (in ar- 
bitrary units) for the double-band case for parameters (a) 

{Sgl = -1, ^32 = 0), (b) (^91 = -2, Sg2 = 0), (C) (Sgl = S, 
Sg2^0). 
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FIG. 5. The absorption and dispersion spectra (in arbi- 
trary units) for the single-band case for parameters 7 = 1, 
Sg = 0. The sohd curve is the absorption profile (— lm[x((5)]) 
while the dashed curve the dispersion profile (Re[x('5)]). 
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FIG. 6. The absorption and dispersion spectra (in ar- 
bitrary units) for the double-band case for parameters (a) 
(&,i = -1, S,2 = 1), (b) (5,1 = -2, 5,2 = 2), (c) (5,1 = -3, 
3,2^3). 
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